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ABSTRACT 


An  experimental  investigation  of  the  restart  process 
of  an  oxygen-Plexiglas  hybrid  rocket  demonstrated  that 
preheating  the  fuel  (as  a  consequence  of  a  previous  igni- 
tion and  a  temporary  extinguishment)  significantly  in- 
creases the  rate  of  chamber  pressurization  and  produces 
regression  rate  overshoots  during  reignition.   Higher 
rates  of  chamber  pressurization  measured  during  the  restart 
transient  imply  faster  instantaneous  regression  rates  of 
the  fuel  during  restart.   Transient  periods  following  the 
initial  ignition  and  the  restart  after  a  two-second  shut- 
down were  observed  for  four  experimental  tests  which  were 
conducted  at  two  oxidizer  flow  rates,  0.078  lbm/sec  and 
0.039  lbm/sec,  and  two  chamber  pressures,  80  psig  and  320 
psig.   Thermocouples  made  from  0.001  inch  diameter  chromel- 
alumel  wire  were  embedded  in  the  fuel  to  record  subsurface 
temperature  histories. 

A  mathematical  model  of  the  thermal  processes  in  the 
fuel  and  the  events  that  occur  during  an  experimental 
firing  was  developed  to  calculate  the  regression  rate  tran- 
sients.  Quantitative  agreement  with  experimental  results 
was  obtained  by  making  nominal  corrections  to  the  calcu- 
lated convective  and  radiative  heating  rates .   A  parametric 
study  investigated  the  influence  of  the  fuel's  thermal 
characteristics  (thermal  conductivity,  surface  temperature, 
and  heat  of  gasification)  and  the  shut-down  duration  on  the 
transient  responses.   The  calculated  results  showed  that 
instantaneous  regression  rate  overshoots  were  as  high  as 
29%  when  high  energy  fuels  are  reignited. 

The  results  of  this  study   indicate  that  the  transient 
behaviors  of  pressure  and  reyression  rate  during  reignition 
after  shut-down  (the  duration  of  which  has  been  varied  over 
a  range  of  0.9  to  4.0  sec)  can  be  anticipated  and  predicted 
once  the  regression  characteristics  of  individual  fuel/oxi- 
dizer  combinations  have  been  established.   Thus,  extensive 
experimental  firings  of  full  scale  rocket  motors  will  not 
be  required  to  adequately  predict  pressure-time  histories 
and  fuel  consumption  following  reignition  if  proper  atten- 
tion is  given  to  understanding  the  thermal  distributions 
during  shut-down  periods . 
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CHAPTER  I 


INTRODUCTION  TO  THE  CONCEPT  OF  THE 
HYBRID  COMBUSTION  MOTOR 


I -A:   Preliminary  Remarks 

Research  and  development  programs  have  advanced  the 
hybrid  rocket  motor  state-of-the-art  to  the  point  that  it  is 
at  times  competitive  and,  in  several  applications,  superior 
to  either  conventional  liquid  or  solid  propellant  rocket 
motors.   The  application  of  hybrid  rocket  motors  in  ad- 
vanced propulsion  systems  is  imminent.   The  hybrid  rocket 
features  the  advantages  of  a  high  specific  impulse  in  a 
storable  system,  a  wide  range  of  thrust  modulation,  stop- 
start  capabilities,  and  the  safety  of  separate  fuel  and 
oxidizer  storage.   A  hybrid  rocket  consists  of  a  gaseous 
or  vaporized  oxidizer  passing  over  a  solid  fuel.   The  high 
specific  impulse  is  obtained  through  the  use  of  special 
high  energy  ingredients  in  the  solid  fuel  and  by  using 
high  energy  fluid  oxidizers.   The  inverse  of  this  con- 
figuration has  been  demonstrated;  the  inverse  hybrid 
consists  of  a  solid  oxidizer  and  liquid  or  gaseous  fuel. 
Conceptual  representations  of  both  the  conventional  and 
inverse  hybrids  are  shown  in  Fig.  1. 

The  combustion  processes  of  a  hybrid  rocket  may  be 
considered  to  consist  of  a  reacting  boundary  layer  with 
heat  and  mass  transfer  as  well  as  chemical  reactions  occuring 
simultaneously.   The  combustion  model  for  an  oxidizer  flowing 
over  a  flat  plate  of  fuel  is  shown  in  Fig.  2.   After  ignition 
is  achieved,  the  injected  gaseous  oxidizer  and  the  fuel  vapor 
(vaporized  from  the  solid  fuel  by  the  heat  feedback  from  the 
flame  zone)  are  mixed  in  the  boundary  layer  and  form  a  tur- 
bulent diffusion  flame  zone. 

The  combustion  rate  of  the  hybrid  rocket  motor  ingre- 
dients depends  on  a  complex  interplay  of  the  following: 
(1)  total  mass  flux  of  the  oxidizer,  (2)  total  mass  flux  of 
the  fuel,  (3)  convective  and  diffusive  mass  transfer  of  the 
reactants  and  products  in  the  turbulent  boundary  layer, 
(4)  convective  and  radiative  heat  transfer  from  the  flame 
to  the  surface  and  (5)  motor  and  chamber  configurations 
that  may  promote  or  retard  the  combustion  processes.   The 
resulting  interactions  that  exist  between  fuel  regression 
rate  and  fuel  geometry,  wherein   the  point  regression  rate 
is  a  function  of  the  upstream  regression  rates  and  the  flow 
conditions,  posed  serious  problems  to  the  development  of 
hybrid  motors.   No  simple  regression  rate  expression,  such 
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as  Tx  elf   as  applied  to  conventional  solid  propellants, 
could  be  prescribed  by  previous  hybrid  researchers;  rather, 
a  complex  combustion  model  is  needed  to  predict  the  rate 
of  regression  of  the  solid  phase  as  a  function  of  the 
geometry  of  the  motor,  the  flow  rate  of  the  oxidizer,  the 
operating  pressure,  and  the  characteristics  of  the  fuel. 
Because  of  the  complexities  of  the  combustion  process, 
considerable  difficulty  is  being  experienced  in  the  devel- 
opment of  rocket  motor  configurations  that  produce  uniform 
regression  of  the  solid  fuel. 

A  large  portion  of  the  hybrid  rocket  research  has 
focused  on  the  following  two  aspects  of  this  surface  regres- 
sion problem:   (a)  a  physical  model  describing  the  reaction 
in  the  boundary  layer  as  a  function  of  the  physical  proper- 
ties of  the  propellant  ingredients  to  yield  a  thermal 
driving  force  and  (b)  a  mathematical  model  of  the  regression 
rate  that  encompasses  the  thermodynamic  properties  of  the 
propellant  and  the  fluid  dynamics  of  the  flow.    The  devel- 
opment of  these  regression  rate  models  has  coincided  with 
advances  made  in  the  fundamental  investigations  of  the 
hybrid  combustion  mechanism.   The  results  of  more  recent 
studies  are  reviewed  in  Chapter  II. 

I-B :   Problems  in  Connection  with  Ignition,  Shut-down  and 
Restart 

While  researchers  have  made  significant  advances  in 
fundamental  investigations  of  the  combustion  mechanism  2,3/ 
the  development  of  more  sophisticated,  high  energy  propel- 
lant combinations  with  additives  and  high  metal  loading  in 
polymer  fuels  has  produced  new  developmental  problems.   Of 
particular  interest  to  practical  hybrid  configurations  are 
the  problems  that  might  occur  during  ignition,  shut-down, 
and  restart  periods.   They  may  arise  as  a  result  of  one 
or  more  of  the  following  factors: 

1.  A  high  loading  of  metal  powder  used  in  some  hy- 
brid fuels  imparts  a  high  thermal  conductivity 

to  the  fuel.   Thus,  more  heat  must  be  transferred 
to  the  subsurface  regions  to  achieve  ignition. 
As  a  result,  pronounced  delays  in  achieving  steady- 
state  operation  are  encountered. 

2.  When  partially  oxidized  fuel  charges  are  employed 
for  enhancement  of  specific  impulse  or  for  enhance- 
ment of  regression  rate,  the  exothermic  reactions 
induced  by  the  burning  process  during  the  normal 
thrust  period  continue  rather  vigorously  after 
oxidizer  shut-off,  and  so  sharp  thrust  termination 
becomes  uncertain. 
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3.  The  slow  regression  rate  characteristic  of  most 
hybrid  fuels  produces  a  thick  heated  layer  that 
remains  on  the  fuel  surface  after  shut-off  of 
the  oxidizer  at  thrust  termination,  and  this 
residual  heat  can  produce  continued  charring, 
surface  reaction,  and/or  hot  fuel  vapor  evolu- 
tion for  some  time  afterward. 

4.  If  the  state  of  the  surface  is  altered  consid- 
erably as  a  result  of  a  continued  reaction  and 
residual  heating  after  combustion  termination, 
the  restart  process  may  be  adversely  affected. 
For  example,  if  a  deep  heated  region  is  pro- 
duced in  the  fuel,  significant  overpressures 
may  result  during  the  reignition  transient. 

A  complete  understanding  of  these  factors  and  the 
ability  to  predict  their  influence  on  the  ignition  and 
restart  processes  is  a  prerequisite  to  the  design  of  a 
high  performance,  rapidly  responding  hybrid  rocket. 
Overpressurization  and  excessive  rates  of  pressurization 
during  the  reignition  transient  period  may  lead  to  in- 
creases in  the  weight  of  structural  components  and 
erratic  thrust  versus  time  histories. 

In  particular,  the  restart  or  reignition  process 
requires  additional  considerations  as  a  result  of  the 
thermochemical  processes  and  thermophysical  changes  that 
may  occur  in  the  fuel  binder  during  the  shut-down  period. 
Because  both  the  total  redistribution  of  heat  (by  con- 
vection, conduction,  and  radiation)  throughout  the  motor 
and  the  degree  of  preheatedness  of  the  fuel  are  functions 
of  time  during  the  shut-down  duration,  a  formidable  prob- 
lem evolves  in  predicting  the  pressure  transients  during 
the  restart  process. 

Out  of  this  description  emerge  the  problems  that 
require  both  theoretical  and  experimental  considerations 
that  are  essential  to  gain  an  understanding  of  the  igni- 
tion and  restart  processes.   Also,  because  of  the  con- 
tinuing trend  toward  larger,  more  sophisticated,  and 
more  costly  hybrid  motors,  experimental  results  that  can 
be  applied  to  practical  hybrid  rockets  must  be  achieved 
using  laboratory  scale  rocket  motors.   Therefore,  a  study 
has  been  conducted  to  determine  the  effects  of  the  fol- 
lowing factors  on  the  restart  process:   (1)  temperature 
profiles  in  the  fuel,  (2)  heat  redistribution  to  the 
fuel  from  the  motor  following  oxidizer  termination  and 
(3)  the  duration  of  the  shut-down  interval. 
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CHAPTER  II 


LITERATURE  SURVEY 


No  simple  regression  rate  law  is  available  to  the  hy- 
brid rocket  designer  and  analyst.   Instead,  previous  re- 
searchers have  formulated  difficult  to  use,  qualitative 
expressions  that  describe  the  processes  occurring  in  the 
reacting  boundary  layer  adjacent  to  the  regressing  fuel. 
Numerous  studies  in  which  both  simple  and  complex  ap- 
proaches to  the  reacting  boundary  layer  have  provided 
various  regression  rate  expressions.   The  purpose  of  the 
following  discussion  is  not  to  give  a  complete  account  of 
all  the  previous  hybrid  regression  rate  studies  (this  has 
been  expertly  done  in  Ref.  1)  but  rather,  to  indicate  to 
the  reader  the  general  direction  and  scope  of  previous 
research  on  regression  rate  models  for  hybrid  rockets. 
In  Chapter  III,  a  discussion  of  the  factors  influencing 
the  problems  connected  with  shut-down  and  restart  will 
be  treated. 

II-A:   Hybrid  Combustion  Models 

The  nature  of  combustion  within  a  hybrid  rocket, 
wherein  combustion  is  controlled  by  convective  and  radia- 
tive heat  transfer,  rather  than  a  direct  reaction  between 
the  oxidizer  and  the  fuel,  suggests  that  a  combustion 
model  could  be  similar  to  that  of  a  turbulent  diffusion 
flame  within  the  boundary  layer  with  mass  addition.   As 
a  result  of  this  similarity,  the  sought  after  regression 
rate  law  describing  such  a  combustion  model  would  natu- 
rally contain  terms  coupling  it  to  the  aerothermodynamics 
of  the  boundary  layer.   Attention  has  therefore  been  fo- 
cused on  the  chemically  reacting  boundary  layer.   Lees 
presented  a  comprehensive  study  of  convective  heat  trans- 
fer from  reacting  boundary  layers  by  considering  the 
influence  of  fuel  mass  addition  from  the  wall. 

Marxman  and  Gilbert6  applied  Lees'  generalized 
category  of  a  reacting  boundary  layer  to  the  combustion 
process  in  a  hybrid  rocket.   They  found  that  if  Le=l 
when  the  boundary  layer  is  laminar,  and  if  Reynold's 
analogy  is  valid  when  the  boundary  layer  is  turbulent, 
the  regression  rate  is  proportional  to  heat  flux  reaching 
the  surface  from  the  turbulent  boundary  layer  adjacent  to 
the  fuel.   The  result  was  a  regression  rate  equation  based 
on  an  ideal  turbulent  boundary  layer.   The  combustion 
model,  as  shown  in  Fig.  3,  assumes  a  thin  flame  combustion 
layer  above  the  solid  fuel.   Marxman  and  Gilbert's 
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regression  rate  law  illustrates  the  coupling  between  the 
fluid  dynamics  of  the  flow  and  the  aerodynamics  of  the 
boundary  layer.   The  complete  equation  is  given  below:* 

i.KGRja   Cjj   Uj   tth   o-6»(€a~£  -  ec,7w   ) 

Three  important  parameters  are  included  in  the 
equation.   The  first  parameter  is  the  blowing  coefficient 
ratio,  C«/cHo    ,  relating  the  coupling  between  regression 
rate  and  boundary  layer  flow  by  accounting  for  the  partial 
blockage  of  heat  transfer  to  the  fuel  surface  as  a  result 
of  mass  addition  from  the  fuel  surface.   The  second  pa- 
rameter is  the  independent  variable,  aH/&  H   ,  coupling  the 
regression  rate  equation  to  a  thermodynamic  parameter 
dependent  on  only  the  fuel-oxidizer  combination  and  its 
gasification  process.   The  third  parameter  is  the  velocity 
ratio,  Ue/Ug   ,  leading  to  a  solution  locating  the  flame 
position  relative  to  the  boundary  layer  thickness.   Marxman 
and  Gilbert  have  also  included  a  radiation  term  to  include 
any  additional  heat  transfer  as  a  result  of  radiative 
energy.   However,  during  normal  operating  conditions,  the 
radiative  addition  is  likely  to  be  small  in  comparison  to 
the  convective  heat  transfer.   Applying  Eq.  1,  less  the 
radiation  term,  to  an  oxygen-Plexiglas  system,  Marxman 
and  Gilbert  have  correlated  regression  rates  over  a  wide 
range  of  oxidizer  mass  fluxes  and  Reynolds  numbers  and 
obtained  agreement  with  experimental  values  within  a 
factor  of  two.   However,  as  noted  by  Marxman  and  Gilbert, 
experimental  and  theoretical  work  is  still  needed  to 
clarify  the  influences  of  the  parameters  in  the  equation. 

It  is  convenient  at  this  time  to  introduce  and  define 
the  mass  transfer  number  &?(Ue/</g)  (aH/aM)    •   This 
parameter  plays  a  dual  role  in  hybrid  combustion.   Initially, 

fi  may  be  considered  as  a  thermodynamic  constant  charac- 
terized by  the  physical  properties  of  the  solid  and  the 
gas  phase  of  the  particular  propellants.   Further,  it  has 
also  been  shown  that  B=  Pr    wfl  is  the  similarity  parameter 
of  a  boundary  layer  with  mass  injection.   Thus,  if  ft  is 
constant  along  the  wall,  the  velocity,  species  concen- 
tration, and  enthalpy  profiles  are  similar  in  the  boundary 
layer.   Since  the  Prandtl  number  of  the  gas  is  near  unity, 

B'  and  B  ar^  not  only  similar,  but  also  approximately  equal. 


*A11  symbols  are  defined  in  the  nomenclature. 
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Thus,  O  serves  to  couple  the  regression  rate  equation  to 
a  particular  propellant  combination  and  the  fluid  dynamics 
of  the  flow.1 

An  empirical  analysis  of  the  regression  rate  equation 
was  performed  by  Anderson.7   Anderson7  applied  the  Dittus- 
Boelter  heat  transfer  equation  to  turbulent  flow  over  a 
fuel  sample  in  a  blast  tube.   Using  the  Dittus-Boelter 
equation 

A  (2) 

and  a  constant  Prandtl  number  of  unity 


X 


the  heat  transfer  coefficient  becomes 

0.8 


(3) 
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Anderson  then  applies  the  ideal  gas  law  to  express  f   ,  "\T. 
and//  (and  thus  ReJ  in  terms  of  temperature,  T,  pressure, 
P,  and  Mach  Number,  M.   The  heat  transfer  rate  is  equal 
to  the  convective  heat  transfer  coefficient  multiplied  by 
the  temperature  difference,  Tf  -  Tw. 

o.ft 


The  regression  rate  will  be  directly  proportional  to  the 
amount  of  heat  reaching  the  surface  and  inversely  pro- 
portional to  the  heat  required  to  vaporize  the  solid. 
Hence,  a  simplified  regression  rate  equation  is  given: 


(5) 


i 


aH 


0.8 
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In  summary,  Anderson  related  the  regression  rate 
equation  to  Mach  number,  pressure,  and  the  heat  of  abla- 
tion of  the  solid  fuel.   The  equation  is  independent  of 
the  oxidizer  choice,  and  his  results  for  the  regression 
rate  of  Plexiglas  as  a  function  of  the  Mach  number- 
pressure  produce,  M'P,  is  shown  in  Fig.  4.   Anderson 
believes  that  this  practical  expression,  in  terms  of 
easily  obtained  parameters,  is  a  good  approximation  over 
a  wide  range  of  the  oxidizers  and  is  therefore  a  useful 
engineering  tool  for  the  design  of  hybrid  motors. 

In  an  extension  of  Anderson's  work,  Hurt  and 
Anderson**  arrived  at  Eq.  6  by  an  approach  that  included 
experimental  and  theoretical  examination  of  the  connec- 
tive heat  transfer  through  the  laminar  sublayer  protion 
of  the  turbulent  boundary  layer.   The  laminar  sublayer 
normally  offers  the  greatest  resistance  to  heat  flow. 
Hurt  and  Anderson^  postulated  that  the  ablation  of  the 
solid  fuel  is  controlled  by  the  thickness  of  the  viscous 
sublayer  and  that  a  heat  transfer  coefficient  with  zero 
fuel  mass  addition  could  be  chosen  and  multiplied  by  a 
correction  factor  to  account  for  the  effect  of  fuel  mass 
addition  at  the  wall. 

The  pressure  dependence  of  both  metalized  and  non- 
metalized  hybrid  fuels  has  received  considerable  atten- 
tion. 9, 10, 11   Smoot  and  Price  concluded^ ,10  that  in 
regions  of  lower  oxidizer  flow  rates  the  regression  rates 
were  independent  of  pressure  and  increased  as  the  0.8 
power  of  specific  total  flow  rate.   However,  at  high  flow 
rates,  the  regression  rate  was  nearly  independent  of  flow, 
but  increased  markedly  with  pressure.   Thus,  their  expperi- 
mental  results  are  in  good  agreement  with  those  predicted 
by  the  turbulent  heat  transfer  model  for  only  the  lower 
flow  rates.   It  was  therefore  postulated  by  Smoot  and 
Price^'lO  that  rate  limiting  chemical  kinetic  processes 
were  the  cause  of  the  observed  pressure  dependence. 

It  is  important  to  mention  that  the  approaches  to 
modeling  a  regression  rate  law,  whether  it  is  the  turbu- 
lent boundary  layer  combustion  model  or  the  M»P  correlation, 
have  primarily  been  directed  at  the  steady-state  condi- 
tion.  That  is,  the  regression  rate  is  directly  proportional 
to  the  heat  reaching  the  surface.   Therefore,  neither 
approach  has  been  extended  to  consider  a  transient  period 
when  the  instantaneous  regression  rate  is  changing  to 
approach  its  steady-state  value.   A  transient  regression 
rate  may  arise  as  a  result  of  any  of  the  following:  (1) 
after  ignition,  the  burning  conditions  approach  the  steady- 
state  value,  (2)  a  throttling  maneuver  when  a  rapid  change 
in  oxidizer  flow  is  followed  by  an  increasing  or  decreasing 
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regression  rate  and  (3)  termination  of  the  oxidizer  flow 
when  the  regression  rate  approaches  zero. 

II-B:   Time  Dependent  Operation 

An  analytical  investigation  of  the  time  dependent 
transient  period  has  been  performed  by  Marxman,  Wooldridge 
and  Muzzy.  2   Their  results  are  directed  at  only  the  tran- 
sient period  following  initiation  of  heat  to  the  surface 
(ignition) .   This  transient  period  is  related  to  the  time 
required  to  achieve  a  steady-state  thermal  profile  in  the 
fuel.   During  this  period,  the  regression  rate  is  not 
proportional  to  the  heat  flux  received  at  the  surface  since 
the  thermal  profile  is  developing  to  the  steady-state  con- 
ditions.  During  ignition  transients  the  ratio  of  the  heat 
conducted  into  the  solid  to  the  heat  that  gasifies  the  fuel 
is  higher  than  for  the  steady-state  conditions.   As  a 
result,  regression  rate  equations  based  on  steady-state 
thermal  balances  at  the  gas-solid  interface  do  not  apply 
during  transient  conditions.   Thus,  it  is  necessary  to 
account  for  the  developing  thermal  profiles  in  the  fuel 
in  order  to  predict  the  instantaneous  regression  rate. 

Marxman,  et  al .  continued  their  analysis  to  predict 
the  depth  of  the  thermal  profile  as  a  function  of  time. 
The  steady-state  penetration  depth,  Sao      t    which  the  pro- 
tile  will  reach  is  proportional  to  the  thermal  diffusivity 
and  inversely  proportional  to  the  steady-state  regression 
rate  as  approximated  by  Eq.  7. 


A.=3  * 


Poo  (7) 


The  nomenclature  associated  with  the  steady-state  thermal 
profile  is  shown  in  Fig.  5. 

As  a  result  of  their  analysis  the  transient  behavior 
of  the  regression  rate  can  be  expressed^s  a  function  of 
the  parameter  C^  (Ts-T0  )/*  H    (assuming  r^/K  is  constant) 
as  shown  in  Fig.  6.   Consequently,  the  transient  periods 
are  greater  for  fuels  with  a  low  An  or  a  high  vaporization 
temperature.   However,  the  possibility  of  decreasing  the 
transient  period  arises  as  a  result  of  increasing   T  , 
the  initial  fuel  temperature.   This  is  not  considered  by 
Marxman,  Wooldridge  and  Muzzy.   Although   TQ  is  considered 
the  constant  temperature  deep  in  the  fuel,   TQ  may  be 
increased  as  a  result  of  a  previous  firing. 

In  conclusion  then,  the  relationship  established 
between  the  transient  regression  rate  and  the  time  to 
establish  a  thermal  profile  can  be  extended  one  step  further 
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As  the  steady-state  thermal  profile  is  established  and 
more  fuel  mass  is  added  to  the  combustion  process,  the 
rate  of  pressurization  in  the  chamber  increases  and  also 
depends,  therefore,  on  the  time  to  establish  the  steady- 
state  thermal  profile.   This  interdependence  between 
the  temperature  profile  of  the  fuel,  regression  rate  tran- 
sients, and  the  chamber  pressurization  rate  will  be  dis- 
cussed in  more  detail  in  the  next  chapter. 
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CHAPTER  III 

THE  EFFECTS  OF  THE  THERMAL  ASPECTS  OF  A  HYBRID  FUEL  ON 
THE  SHUT-DOWN  DURATION  AND  RESTART  PROCESSES 

III-A:   Influence  of  a  Preheated  Fuel 

The  purpose  of  this  discussion  is  to  examine  the  prob- 
lems that  may  arise  in  predicting  the  re  ignition  process  as 
a  result  of  a  nonuniform,  nonsteady  thermal  profile  in  the 
fuel.   A  search  of  previous  efforts  revealed  that  insuffi- 
cient emphasis  has  been  placed  on  the  factors  influencing 
the  re  ignition  process.1-3   Therefore,  an  experimental  pro- 
gram was  undertaken  to  determine  the  influence  of  a  pre- 
heated fuel  on  the  reignition  process  by  comparing  the 
pressure-time  histories  and  subsurface  temperature  measure- 
ments of  the  ignition  and  reignition  transients  leading  to 
steady-state  operation. 

As  mentioned  previously,  problems  connected  with  shut- 
down and  reignition  may  arise  as  a  result  of  many  factors. 
By  way  of  review,  a  few  of  the  factors  are:   (1)  charring 
of  the  fuel  as  a  result  of  residual  heat  in  the  fuel,  (2) 
continued  exothermic  surface  reactions  after  termination  of 
oxidizer  flow  and  (3)  heated  subsurface  layers  causing  sig- 
nificant overpressur izations  during  the  reignition  tran- 
sients.  Furthermore,  the  ignition  criteria  applied  to  the 
initial  ignition  model  may  become  invalid  for  reignition  as 
a  result  of  nonuniform  heat  transfer  to  the  fuel  during 
shut-down. 

Therefore,  for  the  design  of  advanced  maneuvering  miss- 
iles to  include  accurate,  predictable,  stop-start  operation, 
emphasis  must  be  placed  on  the  thermal  aspects  influencing 
the  shut-down  period  and  the  reignition  process. 

III-B:   Thermal  Aspects  of  the  Shut-down  Period 

The  following  sequence  of  events  provides  a  model  for 
considering  the  thermal  aspects  that  influence  the  shut- 
down period: 

1.  A  hybrid  rocket  motor  is  operating  at  steady- 
state  conditions. 

2.  The  flow  of  oxidizer  into  the  motor  ports  is 
terminated,  causing  chamber  pressure  to  de- 
crease.  As  a  result  of  this,  heat  transfer 
to  the  propellant  surface  is  reduced.   The 
net  effect  is  extinguishment. 
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3.  Heat  contained  in  the  fuel  at  the  surface  is 
conducted  into  the  interior  of  the  fuel. 

4.  The  fuel  may  continue  to  receive  conductive 
heat  energy  from  continued  exothermic  reactions 
in  a  metal  loaded  fuel  even  after  the  cessation 
of  oxidizer  flow.   Also,  the  flow  of  the  hot 
premixed  combustion  gases  remaining  in  the 
chamber  will  contribute  to  the  convective  heat 
transfer  to  the  fuel. 

5.  The  fuel  may  also  receive  heat  in  the  form  of 
radiant  energy. 

Since  the  propellant  surface  is  cool  in  comparison 
to  motor  components,  such  as  the  nozzle,  the  redistri- 
bution of  the  total  heat  energy  will  proceed  toward  an 
equilibrium  value.   However,  this  redistribution  of 
heat  is  both  nonuniform  and  nonsteady.   A  model  for  the 
nonuniform  distribution  of  heat  transfer  to  the  fuel  is 
shown  in  Fig.  7.   As  a  result  of  this  redistribution,  if 
excessive  heat  is  conducted  into  the  subsurface  layers, 
a  thick  layer  of  the  fuel  is  effected.   For  example,  in 
a  70-percent  lithium-fluorine  hybrid  system,  melt  layer 
depths  of  greater  than  0.125  inches  have  been  measured. 14 
Reignition  of  this  type  of  surface  allows  for  much  of 
the  fuel  to  be  stripped  away  from  the  surface  and  not 
vaporized  before  it  is  ejected  from  the  motor.   This 
can  greatly  reduce  combustion  efficiency  and  may  hinder 
or  even  prevent  reignition. 

In  conclusion,  the  degree  of  heating  of  the  fuel  at 
a  given  location  is  dependent  upon  motor  geometry,  tem- 
perature distribution  throughout  the  motor,  characteristics 
of  the  propellant,  and  duration  of  the  shut-down  period. 

III-C:   Thermal  Aspects  to  be  Considered  in  Formulating  A 

Reignition  Model 

There  is  presently  no  model  that  can  be  used  to 
predict  the  prerequisite  igniter  energy  or  delivery  rate 
for  the  reignition  of  a  preheated  hybrid  fuel.  The  purpose 
of  the  following  discussion  is  not  to  develop  such  a  model, 
but  rather  to  incorporate  the  results  gained  from  a  thermal 
study  of  the  shut-down  period  into  a  descriptive  model  that 
can  be  used  to  predict  a  satisfactory  reignition  process  in 
practical  hybrid  rockets. 

As  the  combined  convective  and  radiant  heat  transfer 
to  the  fuel  during  shut-down  is  conducted  into  the  fuel 
and  away  from  the  surface  (as  described  in  the  previous 
section),  a  thermal  wave  travels  through  the  fuel.  Since 
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the  fuel  is  no  longer  burning,  i.e.  surface  is  no  longer 
regressing,  the  heating  penetrates  deeper  into  the  fuel 
than  for  steady-state  conditions. 

The  temperatures  in  the  layers  near  the  surface  will 
decrease  as  the  deeper  sublayers  are  preheated.   An  ex- 
ample of  the  thermal  profile  established  during  the  shut- 
down interval  is  compared  to  a  steady-state  profile  in 
Fig.  8,  curves  1  and  2  respectively.   The  slope  and  depth 
of  curve  2  is  particularly  significant  in  the  reignition 
model  as  they  determine  the  amount  of  igniter  energy  needed 
to  reignite  the  motor  and  also  determine  the  transient 
regression  rate  behavior  prior  to  the  steady-state  con- 
ditions. 

It  seems  evident  that  if  the  same  igniter  energy 
as  applied  on  the  initial  ignition  is  applied  to  the 
preheated  fuel,  the  following  sequence  should  occur: 

1.  Ignition  delay  will  decrease  because  less  energy 
is  absorbed  by  the  fuel  before  it  reaches  igni- 
tion conditions.   The  time  between  the  onset  of 
the  ignition  stimulus  and  the  actual  reignition 
event  must  be  predictable  to  control  accurate 
thrusting  maneuvers. 

2.  As  a  result  of  the  already  preheated  sublayers, 

a  thermal  profile  (curve  3)  is  established  which 
penetrates  deeper  into  the  fuel  than  the  steady- 
state  profile  (curve  2) .   The  preheated  fuel 
requires  less  heat  flux  to  the  surface  to  bring 
the  fuel  to  its  vaporization  point.   Thus  a 
greater  fraction  of  the  heat  flux  can  be  used 
for  vaporizing  the  fuel.   As  a  result,  instan- 
taneous regression  rates  will  be  higher.   The 
transient  behavior  of  the  regression  rate  follow- 
ing reignition  will  rise  more  rapidly.   If  a 
sufficiently  deep  layer  of  fuel  has  been  pre- 
heated, the  regression  rate  may  overshoot  its 
steady-state  value.   Due  to  the  nonuniformity 
of  the  thermal  profile  along  the  fuel,  the 
regression  rates  along  the  fuel  grain  will  not 
be  uniform.   The  overall  effect  is  to  complicate 
the  axial  dependence  of  regression  rate  and  to 
enhance  the  possibility  of  a  local  burn-out  of 
the  fuel,  thereby  prematurely  exposing  the  motor 
case  to  the  hot  combustion  gases. 

3.  If  higher  instantaneous  regression  rates  are 
encountered,  higher  chamber  pressures  following 
the  reignition  stimulus  will  occur.   These  will 
result  from  higher  rates  of  fuel  mass  addition 
to  the  combustion  process.   The  reignition  model 
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must,  therefore,  consider  the  degree  of  preheat- 
edness  of  the  fuel  to  account  for  the  increased 
fuel  mass  addition  in  order  to  anticipate  the 
chamber  pressure  within  structural  limits  of 
the  motor. 

In  summary,  if  the  same  ignition  stimulus  is  applied 
to  a  preheated  fuel,  the  degree  of  preheatedness  influences 
the  reignition  process  as  follows:   ignition  will  occur 
sooner,  instantaneous  regression  rates  will  be  higher,  and 
chamber  pressures  will  be  greater.   An  experimental  program 
must  first  determine  to  what  extent  a  preheated  fuel  in- 
fluences the  restart  of  a  hybrid  motor  and  then  investigate 
the  effects  of  propellant  properties,  shut-down  duration 
and  igniter  stimulus.   An  oxygen-Plexiglas  hybrid  motor 
was  designed  for  this  program  and  is  described  in  Chapter 
IV. 
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CHAPTER  IV 


EXPERIMENTAL  APPARATUS 


IV-A:   Overall  Gas  System 

A  hybrid  rocket  motor  was  designed  and  constructed  for 
experimental  investigations  of  the  ignition  and  restart  pro- 
cesses.  Figs.  9  and  10  show  schematic  drawings  of  the  hy- 
brid motor  and  overall  gas  system  respectively.   The  system 
consists  of  oxygen,  methane  and  nitrogen  sources  and  their 
associated  regulators  and  gauges.   Oxygen  is  stored  in  type 
1A  bottles  containing  244  cu.  ft.  at  2200  psig,  and  its  flow 
is  regulated  from  the  test  cell  by  a  nitrogen  dome  regulator. 
Immediately  upstream  of  the  solenoid  and  flow  metering  ori- 
fice, an  oxygen  reservoir  holds  15  in.   at  the  system  stag- 
nation pressure.   The  pressure  in  this  reservoir  is  regis- 
tered on  the  oxygen  reservoir  gauge  (0-2000  psig)  in  the 
test  cell.   The  oxygen  reservoir  pressure  and  the  orifice 
size  are  matched  to  achieve  the  desired  oxidizer  mass  flows. 

Methane  is  stored  in  type  1A  bottles  at  2200  psig  and 
fed  through  a  two  stage  regulator  to  the  solenoid  in  the 
methane  line.   The  pressure  set  on  the  second  stage  of  the 
regulator  and  the  orifice  downstream  of  the  solenoid  were 
matched  to  achieve  stoichiometry  between  the  ignition  gases, 
methane  and  oxygen. 

The  orifices  used  as  flow  metering  devices  were  cali- 
brated using  a  gas-water  displacement  apparatus.   All  the 
orifices  tested  were  found  to  be  ±3%  of  theoretically  pre- 
dicted flows. 

Nitrogen,  stored  at  375  psig,  was  used  to  purge  any 
residual  burning  and  remove  combustion  gases  from  the  sys- 
tem.  The  motor  is  shown  in  an  exploded  view  in  Fig.  11  with 
its  components  labeled. 

IV-B:   Description  of  the  Hybrid  Motor 

The  fuel  consisted  of  a  hollow  cylinder  of  Plexiglas 
10.0  inches  in  length  and  3.0  inches  in  outer  diameter.   The 
initial  port  diameter  was  1.0  inch  for  all  experimental  fir- 
ings.  Final  port  diameters  were  always  less  than  1.2  inches. 
Gaseous  oxygen  was  used  as  the  oxidizer.   Experiments  were 
conducted  at  two  oxygen  flow  rates,  0.0785  lbm/sec  and  0.0392 
lbm/sec.   Experiments  were  conducted  over  a  range  of  steady- 
state  chamber  pressures  between  6Q  and  700  psig.   Ignition 
was  achieved  by  a  methane-oxygen  torch  igniter  mounted  on 
the  forward  housing.   Igniter  durations  were  varied  between 
50  msec  and  500  msec.   A  plenum,  or  mixing  chamber,  located 
aft  of  the  grain,  was  designed  to  receive  the  hot  combustion 
products  and  eject  them  out  the  nozzle.   The  motor  in  Fig. 12 
is  shown  assembled  and  mounted  for  firing  in  the  test  cell. 
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A  normal  test  firing  consisted  of  one  or  two  5  to  10  sec 
successive  firings  with  a  predetermined  shut-down  duration. 
A  photograph  of  the  hybrid  rocket  in  steady-state  operation 
is  shown  in  Figure  13. 

IV-C:   Timing  System 

Four  automatically  resetting  timers  were  integrated  on 
the  firing  console.   Two,  0-30  sec,  ATC  timers  (series  305D) 
provided  selection  of  the  firing  duration  and  shut-down 
duration.   The  hybrid  firing  duration  timer  controlled  power 
to  the  solenoid  in  the  oxygen  line.   At  the  end  of  the  fir- 
ing duration,  the  timer  would  signal  the  solenoid  to  close 
and  initiate  the  shut-down  timer.   The  firing  duration  timer 
would  reset  automatically  for  a  duplicate  second  firing  dur- 
ation. 

Two  other  timers  were  available  for  setting  the  ig- 
niter durations.   For  single  firing  runs  the  #2  cam  of  the 
firing  sequencer  could  provide  the  desired  igniter  dura- 
tion.  However,  accurate  positioning  of  the  cam  proved  to  be 
time  consuming  and  not  reproducible.   For  two  successive 
firings,  an  initial  igniter  duration  could  be  set  on  the 
sequence  cam,  while  the  restart  igniter  duration  could  be 
set  on  the  igniter  duration  timer  mounted  on  the  console. 
This  timer  had  a  0-1000  msec  range.   Single  firing  igniter 
durations  sould  also  be  set  on  the  igniter  duration  timer. 
The  timer  controlled  the  opening  and  closing  of  the  sole- 
noid valve  in  the  methane  line.   Figure  14  shows  the  electri- 
cal schematic  drawing  for  the  entire  system. 

IV-D:   Hybrid  Ignition  System 

The  ignition  of  polymethylmethacrylate  (commonly  re- 
ferred to  as  Plexiglas  or  PMM)  was  achieved  with  a  methane- 
oxygen  gas  torch.   The  igniter  system  allowed  both  igniter 
mass  flow  rates  and  igniter  durations  to  vary  as  a  controlled 
experimental  parameter.   Methane  was  metered  through  a  flow 
metering  orifice  in  approximate  stoichiometric  proportion  to 
the  experimental  oxygen  flow  rate.   The  igniter  gases  were 
then  mixed  and  ignited  in  the  forward  housing  and  then  in- 
jected into  the  motor  port. 

The  actual  ignition  sequence  occurred  as  follows:   (1) 
two  spark  plugs  were  actuated  for  500-700  msec,  (2)  methane 
valve  was  opened,  (3)  oxygen  valve  was  opened,  (4)  hot  ignit- 
er gases  entered  the  motor,  (5)  ignition  was  achieved, 
(6)  methane  flow  was  terminated,  and  (7)  constant  oxygen 
flow  was  supplied  for  the  duration  of  the  firing. 

IV-E:   Description  of  the  Hybrid  Fuel 

Plexiglas  has  been  used  extensively  by  other  research- 
ers,  5'   "'   ''    '     before  experimenting  vith  high  energy 
fuels.   As  a  fuel,  PMM  provides  a  safe,  practical,  inexpensive 
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laboratory  model.   Results  obtained  from  this  model  can  be 
used  as  a  guide  for  the  development  of  analytical  models  for 
the  prediction  of  ignition  and  restart  transients  of  other 
hybrid  engines. 

A  cylindrical  shape  was  chosen  for  the  motor  configura- 
tion, rather  than  a  normal  flat  slab  model,  because  a  cylin- 
drical grain  provides  a  realistic  simulation  to  practical 
hybrid  combustion  motors.   However,  there  are  definite  dis- 
advantages inherent  to  the  cylindrical  configuration.   They 
are:   (1)  photography  of  the  burning  mechanism  is  hampered 
by  focusing  through  a  burning  surface  and  (2)  measured  re- 
gression rates  on  curved  surfaces  are  inaccurate.* 

IV-F:   Components  of  the  Aft  Housing 

The  aft  housing  receives  the  hot  combustion  products 
from  the  oxidizer  and  Plexiglas  reaction.   In  order  to  pre- 
vent high  heat  transfer  rates  to  the  pressure  transducer 
housing,  a  high  density  graphite  insert  was  inserted  in  the 
mixing  chamber.   The  graphite  insert  was  made  in  three  sec- 
tions.  The  first  section,  adjacent  to  the  grain,  had  a 
slight  taper  to  provide  a  gradual  change  from  the  port  diam- 
eter to  the  mixing  chamber  diameter.   Without  this  insert, 
complex  flow  patterns  and  undesired  heating  eddies  would  oc- 
cur immediately  aft  of  the  grain. 

Motor  nozzles  were  fabricated  from  copper  to  prevent 
any  measurable  erosion  in  the  throat  section.   Nozzles  were 
designed  with  a  slight  convergent  section  leading  to  a  pre- 
cise throat  diameter.   Required  pressure  seals  were  main- 
tained with  rubber  O-rings. 

IV-G:   Instrumentat  ion 

Chamber  pressure  measurements  were  recorded  during  all 
experimental  firings.   A  Dynisco  pressure  transducer,  Model 
PT-76,  was  mounted  over  the  mixing  chamber  forward  of  the 
nozzle  blocks.   A  1/4  inch  thick  disk  of  General  Electric 
RTV  580  rubber  was  bonded  to  the  diaphram  to  protect  the 
sensitive  diaphragm  from  contact  with  the  hot  combustion 
gases . 

The  same  model  transducers  were  also  placed  downstream 
of  the  solenoid  valves  in  the  oxygen  and  methane  flow  lines 
for  initial  tests.   These  pressure  recordings  were  used  to 
determine  flow  rates  and  opening/closing  times  of  the  sole- 
noid valves.   The  oxygen  solenoid  valve  fully  opened  and 
closed  in  2  5  msec.   The  methane  solenoid  valve  opened  to  the 
full  position  in  40  msec  and  closed  in  2  5  msec.   The  times 
measured  for  firing  duration  and  igniter  duration  were 
measured  from  the  time  the  oxygen  and  methane  valve  received 


*  The  total  regression,  AR  ,  was  measured  with  an  accuracy 
of  i-0.0005  inches  which  gave  reliable  regression  rate  meas- 
urements to  -1%. 
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their  respective  electrical  signals  to  open.   The  signals 

to  open  and  close  the  valves  were  recorded  by  the  visicorder. 

The  RTV  shielded  pressure  transducer  had  a  natural  fre- 
quency of  22000  cps.  2 5^-  was  reliable  to  a  frequency-amplitude 
response  to  8000  cps.     The  transducers  were  driven  with 
B  &  F  Transducer  Conditioning  Modules.   The  pressure  trans- 
ducers were  initially  calibrated  using  the  Amthor  0-2000 
type  472  Deadweight  Tester.   Each  experimental  firing  was 
accompanied  by  a  pressure  calibration  of  the  motor  with  a 
regulated  nitrogen  supply. 

Some  experimental  firings  were  instrumented  with  fine 
bare  wire  thermocouples  embedded  in  the  fuel.   These  thermo- 
couples were  made  from  chromel-alumel  wire  0.001  inches  in 
diameter.   A  detailed  explanation  of  the  technique  developed 
to  embed  these  thermocouples  will  be  given  in  Appendix  I. 

Both  pressure  signals  and  thermocouple  outputs  were  fed 
to  Dana  Model  3520  D.  C.  Amplifiers.   At  the  Honeywell  Model 
1508  Visicorder,  these  signals  were  matched  with  Honeywell 
galvonometers  (Model  1650)  and  were  recorded  at  a  visicorder 
speed  of  12  in/sec  with  timing  marks  every  10  msec. 

Thermocouple  outputs  were  calibrated  prior  to  the  fir- 
ing by  using  a  potentiometer  to  input  a  desired  voltage  to 
the  galvonometers.   Amplifiers  could  then  be  adjusted  to 
match  the  desired  deflection  per  millivolt.   Most  thermo- 
couple outputs  were  calibrated  for  1/2  inch  deflection  per 
3  millivolts. 

High  speed  photography,  using  a  16mm  Wollensak  Fastex 
camera,  provided  observations  of  the  ignition  and  restart 
phenomena.   High  speed  films  can  locate  the  point  of  igni- 
tion, direction  of  flame  spreading,  and  the  rate  of  flame 
propagation  over  the  surface. 

Films  taken  at  1500  and  2000  frames  per  second  were 
electronically  sequenced  to  be  at  operational  speeds  at  the 
time  of  ignition.   Two  neon  marking  lights  were  used  to 
coordinate  recorded  data  and  filmed  observations.   One  mark 
recorded  simultaneously  on  both  the  film  and  visicorder 
served  as  a  reference  mark.   The  second  mark  was  an  oscil- 
lating timing  light  set  at  100  H^  • 
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CHAPTER  V 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 


V-A:   Thermocouple  Measurements 

Thermocouples  were  embedded  in  the  Plexiglas  fuel 
for  five  experimental  firings.   At  each  experimental 
oxidizer  flow  rate,  0.078  lbm/sec  and  0.039  lbm/sec ,  we 
attempted  to  investigate  two  operating  pressures  of 
approximately  80  and  360  psig.   Four  firings,  Runs  17, 
18,  22,  and  23,  were  instrumented  with  four  thermocouples 
per  firing,  each  at  a  different  depth  in  the  fuel  (see 
Table  I).   A  fifth  firing,  Run  19,  was  conducted  with 
five  thermocouples  spaced  axially  along  the  grain  all  at 
the  same  subsurface  depth.   All  five  firings  were  con- 
ducted with  the  same  firing  sequence,  trun=5.0  sec, 
tsj1(^=2.0  sec,  and  tre^=5.0  sec. 

The  temperatures  recorded  during  Runs  17,  18,  22, 
and  23,  were  plotted  versus  time,  in  Figs.  15-18.   Fo- 
cusing our  attention  on  the  shut-down  duration  (5-7  sec) , 
it  is  evident  that  the  temperatures  in  the  subsurface 
layers  continued  to  rise  after  termination  of  the  oxidizer 
flow.   This  temperature  rise  is  caused  by  continued  heat 
conduction  from  higher  temperature  regions  near  the 
surface  and  radiation  from  the  hot  motor  components . 

However,  as  mentioned  previously,  this  redistribu- 
tion of  heat  during  the  shut-down  period  may  not  be 
uniform.   The  results  obtained  from  Run  19,  shown  in 
Fig.  19,  give  supporting  evidence  to  the  nonuniformity  of 
the  heat  distribution.   More  clearly,  Fig.  19  shows  the 
temperatures  recorded  at  1.75,  6.0  and  7.5  inches  from 
the  injector  face  increased  5,  3  and  6  degrees  C  respec- 
tively during  the  shut-down  interval,  while  temperatures 
at  4.5  inches  from  the  injector  face  decreased  3  degrees 
C  during  the  first  second  of  shut-down  and  increased  2 
degrees  C  during  the  last  second  of  shut-down.   In  con- 
clusion, the  results  shown  in  Figs.  15-19  exhibit  the 
nonsteady  and  nonuniform  temperatures  of  the  subsurface 
layers  and  confirm  our  earlier  notions  of  the  nonsteady, 
nonuniform  heat  distribution  during  the  shut-down  interval 

Thermocouple  data  was  also  used  to  measure  the  growth 
of  the  thermal  profile  following  ignition  and  reignition 
stimuli  and  during  the  shut-down  interval.   In  order  to 
experimentally  model  a  thermal  profile,  three  factors 
were  considered:   (1)  instantaneous  regression  rate,ft  t 
(2)  instantaneous  location  of  the  surface  and  (3)  accurate 
knowledge  of  the  thermocouple  location. 
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The  instant  that  the  ablating  surface  passes  each 
thermocouple  bead  can  be  approximated  in  two  ways.   First, 
an  examination  of  voltage  output  from  the  thermocouple 
reveals  a  smooth  voltage  output  while  in  the  solid  region, 
but  as  the  temperature  sensing  bead  enters  the  molten 
liquid  or  gas  phase  region,  the  voltage  output  becomes 
wavy  or  rippled.   Two  examples  taken  from  different 
firings  are  shown  in  Fig.  20.   Second,  plots  of  recorded 
temperature  data  on  a  log  Tx  vs.  time  scale  yield  a  straight 
line  in  the  solid  region  during  steady-state.   However,  a 
noticeable  change  in  the  slope,  normally  in  the  form  of  a 
departure  from  the  straight  line  steady-state  approximation, 
occurs  when  the  thermocouple  records  the  surface  tempera- 
ture, Ts.   Temperature  measurements  and  the  times  the 
thermocouple  beads  reach  the  surface  are  shown  in  Figs. 
21  to  24. 

By  combining  the  two  previous  observations  an  empiri- 
cal regression  rate  can  be  attained  and  applied  to  deter- 
mine the  location  of  the  burning  surface  during  the  entire 
firing.   An  example  of  several  experimental  thermal  profiles 
measured  during  initial  steady-state  operation,  shut-down 
and  reignition  steady-state  operation  are  shown  in  Fig.  25 
for  Run  18.   It  should  be  noted  here  that  the  maximum  error 
in  the  thermocouple  placement  was  0.005  cm  which  represents 
the  sum  of  the  maximum  errors  of  each  individual  measurement 
taken.   For  further  analysis,  the  profiles  will  be  repre- 
sented as  smooth  curves  to  facilitate  any  comparisons. 

The  profiles  shown  in  Fig.  25  particularly  show  the 
overheating  of  the  subsurface  layers  deeper  in  the  fuel 
as  the  shut-down  duration  proceeds.   The  entire  heat  ca- 
pacity of  the  thermal  wave  is  free  to  "cook"  the  deeper 
layers  of  fuel  as  little  or  no  heat  is  used  for  ablation. 
It  also  appears  that  almost  immediately  after  the  reigni- 
tion stimulus  is  applied  at  t  =  7.00  sec,  the  thermal 
profile  established  at  t  =  7.12  sec  is  deeper  than  the 
steady-state  profile  at  t  =  4.   This  would  imply  faster 
instantaneous  regression  rates  because  a  smaller  fraction 
of  the  heat  is  transferred  to  the  subsurface  regions  of 
the  fuel  and  consequently  more  heat  is  available  to  vapor- 
ize more  fuel.   The  results  would  imply  higher  transient 
regression  rates  as  the  preheated  regions  are  consumed. 

Further  comparison  of  the  measured  penetration  depth 
of  the  thermal  profile  during  initial  steady-state  opera- 
tion with  the  steady-state  penetration  depth  predicted  by 
Marxman,  Wooldridge  and  Muzzy-^  in  Eq.  7  show  close  agree- 
ment.  However,  the  measured  thermal  profiles  following 
the  reignition  stimulus  penetrate  deeper  into  the  fuel 
with  temperatures  generally  exceeding  corresponding  steady- 
state  profiles  (Fig.  25). 
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Furthermore,  comparisons  of  the  times  to  reach  steady- 
state  chamber  pressure  (as  tabulated  in  Table  1)  with  the 
times  to  establish  a  steady-state  thermal  profile  agree 
very  closely.   The  time  to  establish  a  steady-state  thermal 
profile  may  be  estimated  graphically  from  a  plot  of  log  T 
vs.  time.   When  the  temperature-time  response  from  a  thermo- 
couple becomes  linear,  a  steady-state  thermal  profile  has 
been  established  as  shown  in  Fig.  26. 

V-B :   Chamber  Pressure  Measurements 

The  major  purpose  of  the  experimental  investigation  was 
to  determine  the  influence  of  a  preheated  thermal  profile  on 
the  reignition  transient  by  comparing  the  pressure-time  his- 
tories of  the  ignition  and  reignition  transients,  respec- 
tively.  These  comparisons  of  pressure-time  histories  were 
made  following  a  two  second  shut-down.   The  pressures  during 
reignition  transients  were  always  greater  than  the  respective 
initial  pressures  during  the  ignition  transient,  as  shown  in 
Figs.  27-30.   However,  both  the  ignition  and  reignition 
processes  approached  similar  steady-state  pressures  in  the 
chamber  as  shown  in  Figs.  31-34. 

Reignition  delays  were  up  to  15%  less  than  ignition 
delays.   Furthermore,  the  times  to  reach  90%  of  steady-state 
pressures  following  reignition  were  one-fourth  of  the  initial 
ignition  times.   Both  ignition  delays  and  the  times  to  reach 
90%  of  steady-state  chamber  pressure,  t     ,  are  tabulated  in 
the  compendium  given  in  Table  1.        .9Pc 

V-C :   Regression  Rate  Measurements 

Average  regression  rates  were  measured  on  each  experi- 
mental firing  by  measuring  the  change  in  web  thickness  per 
firing  duration.   Weight  loss  measurements  were  also  recorded 
and  compared  with  previous  runs  of  the  same  durations . 
Regression  rates  were  measured  at  nine  stations,  space  at 
one  inch  intervals  beginning  from  the  injector  face  of  the 
fuel.   Only  the  regression  rates  measured  at  the  5th  station, 
5  inches  from  the  injector  face,  are  recorded  in  Table  1. 

Experimental  conditions  precluded  the  measurement  of 
instantaneous  regression  rates  during  the  firing  duration. 
From  the  results  of  the  pressure  measurements,  we  can  infer 
that  the  instantaneous  regression  rates  following  reignition 
of  a  preheated  fuel  were  faster  than  their  respective  regres- 
sion rates  following  ignition  of  an  ambient  fuel.   However, 
the  results  obtained  from  measuring  average  regression  rates 
do  provide  meaningful  information   helpful  in  understanding 
other  problems  of  hybrid  combustion. 
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For  example,  total  regression,  AK  ,  rather  than 
regression  rate  X*  ,  was  plotted  versus  the  distance  from 
the  injector  face  for  three  firings  of  the  same  duration, 
but  at  different  operating  pressures,  600,  325  and  80  psig 
respectively.   The  results  in  Fig.  35  show  a  definite  x- 
dependence  on  regression  rate  with  minimum  burning  occur- 
ring at  midstream.   Referring  to  the  same  figure,  (Fig.  35) 
the  pressure  effect  on  regression  rate  is  also  noticeable. 
Runs  10  and  13  at  higher  pressures,  600  and  325  psig  respec- 
tively, have  similar  total  regressions;  Run  15  at  a  lower 
pressure,  76  psig,  has  a  smaller  total  regression.   These 
results  concur  with  the  earlier  mention  of  the  pressure 
dependence  and  the  chemical  kinetics. 

V-D:   Photographic  Observation 

Extensive  photographic  observations  of  the  ignition 
and  reignition  processes  and  the  flame  propagation  during 
the  pressure  transients  were  made.   Although  no  significant 
difference  could  be  detected  between  the  ignition  event  and 
reignition  event,  both  processes  were  characterized  by  long 
flame  spreading  durations  (usually  350  msec)  very  similar 
to  the  long  chamber  filling  intervals.   A  flame  spreading 
interval  during  hybrid  ignition  is  shown  in  Fig.  36. 

The  most  peculiar  result  of  the  photographic  obser- 
vation was  the  unusual  distribution  of  flame  intensity 
throughout  the  motor  core.   The  high  speed  films  show  a 
nonsteady,  almost  pulsating,  nonuniform  flame  distribution 
during  steady-state  operation.   A  sequence  of  photos  in 
Fig.  37  show  the  flame.   Similar  pulsating  flames  were 
reported  by  Courtney^l  in  a  double-slab  02~Plexiglas 
system.   The  nonuniform  intensity  seen  in  both  Figs.  36 
and  37  raises  the  question  of  complete  combustion  along 
the  surface  of  the  fuel.   What  appears  to  be  a  region  of 
"spotty"  higher  intensity  burning  upstream  indicates  that 
upstream  burning  exhausts  the  oxygen  supply.   Accordingly, 
the  oxygen  concentration  in  the  downstream  boundary  layer 
is  not  sufficiently  high  to  burn  the  Plexiglas  at  the 
grain  surface.   One  hypothesis  is  that  hot  combustion 
products  passing  over  the  downstream  regions  of  the 
grain  melts  the  Plexiglas  and  carries  the  melt  to  the 
plenum  chamber  where  it  is  burned. 
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CHAPTER  VI 


A  PARAMETRIC  INVESTIGATION  OF  THE  TRANSIENT  REGRESSION 
RATE  BEHAVIOR  FOLLOWING  A  SHUT-DOWN 


In  order  to  gain  a  better  understanding  of  the  tran- 
sient regression  rate  behavior  on  the  restart  process  of 
a  hybrid  rocket  motor,  a  mathematical  model  was  developed 
to  calculate  the  instantaneous  regression  rates  of  the 
fuel  and  the  temperature  profiles  in  the  fuel.   The  model, 
in  the  form  of  an  extension  and  modification  to  an  existing 
computer  program22,  calculates  the  temperature-time  re- 
sponses at  locations  corresponding  to  the  thermocouples  at 
subsurface  depths.   A  discussion  of  the  computer  solution 
to  the  heat  conduction  equation  used  to  compute  regression 
rates  is  given  in  the  following  section.   Subsequent  to 
that  section,  the  approach  used  in  the  parametric  study 
will  be  treated,  and  the  results  presented. 

VI-A:   Formulation  and  Solution  of  the  Mathematical  Model 

The  mathematical  model  of  the  regressing  fuel  may  be 
considered  in  two  parts:  1)  the  thermal  driving  force  to 
the  fuel  surface  that  gasifies,  decomposes  and  transports 
the  fuel  into  contact  with  the  oxidizing  species  in  the 
chamber  and  2)  the  transient  moving  boundary  heat  conduc- 
tion within  the  fuel. 

VI-A.l:   Heat  Conduction  Equation  for  the  Fuel 

The  equation  governing  transient  heat  conduction  with- 
in a  nonstationary  medium  (in  cylindrical  coordinates)  is 

x     it        an:  :et  3T  (vi_1) 


The  terms  in  Eq.  (VI-1)  have  the  following  physical 
significance : 

I   Time  rate  of  energy  accumulation  and  temperature 
rise  at  a  particular  point 

II   Rate  of  energy  convection  due  to  motion  of  the 
coordinate  system 
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III   Energy  transfer  by  thermal  conduction  in  the  solid 

IV   Term  which  accounts  for  effects  of  varying  cross- 
sectional  geometry  in  the  cylinder  on  energy  trans- 
fer by  thermal  conduction 

V  Rate  of  radiative  heat  transfer  to  the  subsurface 
regions 

The  boundary  condition  at  the  gas-solid  interface  for  the 
fuel  is 


VVV5>«»A(S-)., 


(VI-2 


Where  %~   is  the  convective  heat  flux  from  the  cjombustion 
products  flowing  over  the  fuel.   The  quantity  ^  is  the 
rate  of  heat  release  at  the  gas-solid  interface  due  to 
gasification  of  the  fuel.   For  conventional  fuels  the 
gasification  at  the  gas-solid  interface  is  generally 
considered  to  be  endothermic.   The  quantity^  is  usually 
expressed  in  one  of  two  ways,  in  terms  of  an  Arrhenius 
relation 


v2»"Kia£) 


(VI-3 


or  in  terms  of  a  net  energy  release  per  unit  mass  con- 
verted from  the  condensed  phase  to  the  gas  phase 

V^^  (VI-4 

The  regression  rate  at  the  surface  is  controlled  by 
both  a  linear  pyrolvsis  law, 

VB3e*.p(l^)  (VI_5 

and  a  heat  balance  at  the  surface, 


/A 


(VI-6 
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When  the  pyrolysis  law  (Eq.  VI-5)  describes  the  regres- 
sion rate,  the  temperature  at  the  gas-solid  interface,  T  , 
increases  as  the  heat  flux  to  the  solid-phase  surface  in- 
creases.  As^Tg  increases,  by  Eq.  (VI-5),  a  corresponding 
increase  in  rj>  occurs.   However,  we  chose  to  adopt  a  simpler 
description  of  the  pyrolysis  process,  namely 

Ts  =  Ts'  =  constant  (VI-5)' 

The  value  of  j£  will  now  be  controlled  by  the  other  equations 
of  the  problem,  but  it  will  also  increase  as  the  heat  flux 
increases . 

The  boundary  condition  for  x-**oo  is 

-^—   =  o  (VI_7) 

the  initial  condition  is 

T(ptji0--7i  Cvi-8) 

Because  of  the  low  thermal  conductivity  of  the  fuel  and  the 
rapid  surface  regression,  the  conditions  of  Eq.  (VI-7)  is 
reached  just  0.01  cm  below  the  gas-solid  interface. 

Ignition  is  assumed  to  occur  when  either  the  surface  or 
the  interface  reaches  a  prescribed  ignition  temperature. 

VI-A.2:   The  Thermal  Driving  Force  to  the  Surface  of  the  Fuel 

The  reactions  and  transport  phenomena  within  the  boundary 
layer  adjacent  to  the  fuel  surface,  if  treated  analytically, 
would  be  a  most  formidable  aerothermochemical  exercise.   Sev- 
eral investigators  have  attempted  such  analyses  with  only 
limited  success.^  6,  12  &n   examination  of  the  bulk  processes 
indicate  that  it  may  be  both  attractive  and  meaningful  to  con- 
sider the  convective  heat  transfer  from  the  reacting  gases  to 
the  fuel  using  the  same  functional  relationships  as  those 
developed  for  nonreacting  turbulent  boundary  layers.   There 
are  two  justifications  for  this:   1)  the  work  of  Anderson  and 
Hurt7,  8  demonstrated  that  the  surface  gasification  rates  of 
hybrid  fuels  can  be  correlated  in  terms  of  emperical  turbulent 
boundary  layer  parameters  and  2)  the  analytical  relationships 
for  the  thermal  driving  force  developed  by  Marxman  and  co- 
workers-'-' ^'  l^have  similar  functional  dependencies  on  free 
stream  pressure  and  velocity  as  those  dependencies  on  turbu- 
lent convective  heating. 

Accordingly,  the  thermal  driving  force  used  in  the  model 
was  based  on  the  empirical  relationships^or  high  speed 
turbulent  forced  convection  as  presented  by  Eckert.23 

*   It  should  be  pointed  out  that  we  understood  why  the  simple 
convective  boundary  does  not  fully  account  for  the  diffusion 
blowing,  and  kinetic  processes  within  the  boundary  layer. 
Rather  it  is  used  as  a  systematic  means  of  correlating 
measured  results  within  the  budget  limitations  of  the  study. 
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The  following  general  conditions  were  considered  to 
prevail:   1)  the  boundary  layer  is  steady  and  completely 
turbulent,  2)  an  inviscid  flow  field,  external  to  the 
boundary  layer,  is  imposed  by  the  internal  ballistic 
relationships  of  the  chamber  3)  the  conditions  at  the 
fuel  surface  and  in  the  inviscid  flow  field  are  quasi- 
steady  and  4)  the  heat  conduction  within  the  fuel  is 
unsteady.   In  an  effort  to  account  for  the  effects  of 
temperature  on  the  physical  properties  of  the  gases, 
they  are  evaluated  at  the  following  reference  temper- 
ature 


T  =~T+  +o.s{Tf-T5) 


(VI-9) 


The  convective  heat  transfer  coefficient  based  on  the 
temperature  difference,  HC:  "*  7^   >  is 

h  -  o.ozsg  Q  cp  j>  V  fr      file 

(VI-10) 

The  factor  Cc  is  adapted  to  make  the  calculated  steady- 
state  regression  rates  agree  with  the  measured  steady- 
state  regression  rates. 

VI-A.3:   Solution  to  the  Heat  Conduction  Equation 

The  method   of  explicit  finite  differences  and 
relaxation  were  used  to  solve  the  differential  equations, 
Eq.  (VI-1) .   The  several  terms  in  the  solid-phase  energy 
equation  are  considered  by  writing  the  derivatives  in 
terms  of  the  corresponding  finite  differences.   For 
example, 

3*     ~  At:  -    At  (VI_11) 


and 


VT  _      Tn-f/  -Tn-j 

*X  ~  ZA?C  (VI-12) 


or 


^T  AT" 

»  =  ttl  (vi-13) 
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The  differential  temperature  over  the  time  interval,  At,  is 

(VI-14) 


V   2.  ax   / 


{rC^^^M^^VfW 


Applying  the  techniques  of  explicit  finite  differences, 
the  materials  are  treated  as  a  series  of  elemental  volumes 
through  the  cylinder.   The  mass  of  each  volume  is  consid- 
ered to  be  concentrated  at  a  point  referred  to  as  a  nodal 
point.   Fig.  38  illustrates  an  elemental  volume  and  its 
nodal  point,  n.   The  area  at   x  -Ax/2  is  A^n,  the  area  at 
x  +  &  x/2  is  A out,  and  the  area  at   x   is  A.   The  heat 
transferred  into  the  left  side  of  the  elemental  volume  is 


Qp»  -  K tyi  =A*  JUTo-t^/ax 


(VI-15) 


The  heat  transferred  out  of  the  right  side  is 


Qout * Aout  V** Aut  ^(J^-Tn)/yc 


(VI-16) 


The  temperature   T    is  the  bulk  temperature  of  the  nodal 


volume  whose  heat  capacity  is  £5  pjL   x    (Ain  +  Aout^/2, 
The  heat  flux  at  the  center  of  the  elemental  volume;  i.e., 
at  the  position   x,  is  approximated  as  the  mean  heat  flux 


%M- 


\%IH  +  \o\jt  ) 


(VI-17) 


Substituting  Eqs.  (VI-15)  and  (VI-16)  into  Term  III 
of  Eq.  (VI-14)  yields 


Qc 


ut 


+ 


Q 


K»\      A,MJ«/^ 


At 


(VI-18) 


where 


J 


and 


% 


are  evaluated  at  T 


n 
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Term  IV  in  Eq .  (VI-14)  can  be  expressed  as 

At 


lk- 


J<i 


(VI-19) 


Term  II  accounts  for  motion  with  respect  to  some  fixed 
reference  point.   In  this  analysis  the  solid-phase  material 
is  maintained  at  a  fixed  position  in  space  and  the  burning 
surface  moves  at  the  calculated  burning  rate  relative  to 
that  fixed  position  in  space.   As  the  solid-phase  layers 
are  consumed,  the  nodal  points  within  a  material  are  repos- 
itioned so  that  the  incremental  distances  between  the  nodal 
points,  ^x,   are  equal.   As  the  thickness  of  the  solid 
phase  is  decreased,  the  number  of  nodal  points  in  solid 
phase  is  decreased  so  as  to  maintain  the  stability  of  the 
numerical  solution.   Two  factors  control  the  increment 
sizes  At        and  <£x:   1)  the  desired  accuracy  of  the 
solution;  and,  2)  the  stability  criterion  for  the  finite 
difference  equations.   By  decreasing  the  increment  size 
^x,  the  solution  can  be  made  more  accurate,  but  the  time 
required  for  the  calculations  increases  disproportionately 
since  the  time  increment  At  must  also  be  decreased  to 
maintain  the  stability  of  the  solution.   A  sufficient  con- 
dition for  stability  is  to  insure  that  the  sum  of  all  the 
factors  multiplying  the  temperature  at  a  particular  nodal 
point  is  greater  than  or  equal  to  zero.   In  the  method 
presented  here,  the  time  increment  does  not  vary  as  different 
materials  are  considered.   Thus,  it  is  sufficient  to  con- 
sider the  stability  of  the  critical  nodal  point.   The 
numerical  solution  is  set  up  so  that  the  critical  nodal  point 
for  stability  is  the  nodal  point  at  the  surface  of  the  fuel. 
Accordingly,  stability  is  insured  when 


1  _  z&  - 


(  *■  +  Kw  A*J  -  O  (VI-20) 


Numerous  check  calculations  were  performed  on  all 
portions  of  the  numerical  solution.   In  several  instances 
it  was  possible  to  compare  the  values  calculated  by  the 
numerical  solution  with  the  values  from  closed-form  sol- 
utions and  series  solutions . 

VI-B:   Simulation  of  a  Hybrid  Motor  Test  Firing 

The  purpose  of  the  parametric  study  was  threefold 
(1)  to  compare  the  thermal  profiles  in  the  solid  fuel  estab- 
lished during  steady-state  operations  to  those  established 
during  shut-down  and  reign ition,  (2)  to  compare  computed 
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regression  rate  behaviors  during  restarts  following  dif- 
ferent shut-down  durations  with  the  regression  rate 
behaviors  following  initial  ignition  and  (3)  to  evaluate 
and  compare  the  thermal  response  of  several  hypothetical 
high  energy  fuels. 

The  initial  task  was  to  prescribe  sufficient  input 
conditions  to  the  computer  program  to  simulate  an  experi- 
mental firing.   Run  18  was  chosen  as  the  datum  case  for 
comparison  purposes,  and  the  sequence  that  follows  out- 
lines the  approach  used  to  simulate  the  0  -Plexiglas 
system.  2 

1.  The  physical  properties  of  the  Plexiglas  fuel, 
the  properties  of  the  chamber  gases  and  the 
grain  dimensions  were  input  to  the  program.   No 
parameter  was  input  to  the  initial  datum  case 
to  consider  the  heating  of  subsurface  layers  by 
radiation.   However,  as  further  discussed  in 
Appendix  I,  our  experiments  showed  that  radiant 
heating  to  the  subsurface  layers  of  the  Plexi- 
glas fuel  was  significant  and  must  be  accounted 
for  in  future  investigations. 

2.  The  gas-torch  igniter  was  simulated  by  applying 
a  heat  flux  to  the  surface  of  the  fuel  in  a 
step  function  initiated  at   t  =  0.0.   After 
100  m  sec,  the  step  function  was  followed  by  a 
decreasing  ramp  function  to  zero  at  300  rrrec. 
The  magnitude  of  the  igniter  heat  flux  was  fixed 
at  a  value  that  produced  ignition  delays  that 
agreed  approximately  with  the  experimental  results. 

3.  The  pressure-time  history  recorded  during  Run  18 
was  also  an  input  to  the  program.   From  the 
pressure-time  history,  the  program  calculated  an 
effective  convective  heat  transfer  coefficient, 
h,  and  then  a  heat  flux  to  the  surface,  Q.   The 
heat  flux  to  the  surface  determined  the  regression 
rate.   Small  corrections  were  made  to  computed 
convection  coefficients  to  make  the  regression 
rates  agree  with  experimental  results.   This 
sequence  summarizes  the  inputs  to  the  datum  case; 

a  complete  list  of  input  parameters  is  given  in 
Table  II. 

VI-C:   Calculated  Results  for  a  Plexiglas  Fuel 

The  calculated  results  of  the  datum  case  partially 
simulated  the  experimental  results  of  Run  18.   The  computed 
steady-state  regression  rates  were  within  1%  of  the  measured 
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average  regression  rate.   Ignition  delays  of  0.077  sec 
were  computed  and  were  approximately  equal  to  those  re- 
corded and  photographed  during  the  experiments.   However, 
as  was  indicated  by  a  comparison  of  measured  and  calculated 
thermocouple  responses,  radiant  heating  to  the  subsurface 
regions  should  be  included  in  the  solution  to  the  heat 
conduction  equation. 

An  important  advantage  of  the  computer  simulation  was 
that  thermal  responses  can  be  compared  for  a  range  of  fuel 
properties.   The  thermal  profiles  for  the  nonradiative  datum 
case  are  presented  in  two  forms,  Figs.  39a  and  40a.   In 
Fig.  39a,  the  results  of  the  passage  of  a  thermal  wave 
penetrating  through  the  fuel  during  shut-down  to  preheat 
the  sublayers  above  their  respective  steady-state  tempera- 
tures is  clearly  illustrated.   Curve  1  represents  the 
steady-state  profile  at   t  =  4.0  sec.    Superimposed  on  the 
same  axis,  curve  2  is  the  thermal  profile  2.0  seconds  after 
shut-down  just  prior  to  restart.   The  temperatures  in  the 
layers  adjacent  to  the  surface  have  diminished  slightly; 
however,  the  thermal  wave  has  preheated  the  subsurface 
layers  considerably.   Curve  3  is  the  thermal  profile  at 
t  =  9.0  sec,  less  than  two  seconds  after  reignition.   This 
thermal  profile  is  slowly  approaching  its  original  steady- 
state  condition.   In  Fig.  40a,  the  same  profiles  are  shown 
at  their  instantaneous  surface  locations.   At  t  =  5.0 
sec,  the  reaction  stops  and  the  surface  ceases  to  ablate. 
Therefore,  the  profiles  at  t  =  6.0  and  t  =  7.0  sec  are 
shown  at  the  same  surface. 

The  results  obtained  from  the  temperature-time  his- 
tories at  the  thermocouple  locations  did  not  agree  with 
the  experimental  results  of  Run  18.   The  computed  tempera- 
ture-time histories  lagged  in  response  to  ignition  and  were 
below  the  experimental  values  recorded.   It  was  concluded 
that  the  poor  agreement  resulted  from  considering  only 
convective  heating  in  the  computed  temperatures,  while  the 
experimental  results  were  obtained  by  sensing  convective 
and  radiative  heat  to  the  subsurface  layers. 

Therefore,  two  additional  heating  relationships  were 
considered  to  account  for  the  subsurface  radiant  heating. 
They  were 

(1)   a  radiative  heat  flux  to  the  fuel 
and     (2)   subsurface  heating  by  radiation  using  Beers  Law^ 
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This  introduced  two  additional  parameters  to  the  study, 
the  overall  radiation  interchange  factor,  Af s ,  and  the 
Beers'  Law  parameter,  B. 

With  these  two  appropriate  values  added,  Case  I 
becomes  the  modified  datum  case  considering  a  radiative 
system.   The  temperature-time  responses  of  the  locations 
corresponding  to  the  thermocouples  were  in  reasonable 
agreement.   A  comparison  of  the  measured  and  computed 
temperature-time  responses  for  two  thermocouples  is 
shown  in  Fig.  41.   The  lack  of  better  agreement  between 
the  two  temperature-time  responses  was  probably  due  to  the 
uncertainty  of  the  true  location  of  each  thermocouple  and 
unknowns  associated  with  radiation  effects.   Further 
considerations  of  the  thermocouple  technique  is.  given  in 
Appendix  I. 

A  comparison  of  the  thermal  profiles  for  the  non- 
radiative  datum  case  and  the  radiative  Case  I  is  shown  in 
Figs.  39  and  40.   As  already  described,  Figs.  39a  and  40a 
are  for  the  nonradiative  datum  case,  while  Figs.  39b  and 
40b  show  the  effects  of  radiation  on  the  system.   A 
faster  regression  rate  is  noticeable  by  the  more  advanced 
positions  of  the  instantaneous  surfaces  in  Fig.  40b. 
Also,  the  chemical  reaction  in  the  system  does  not  stop 
at  t  =  5.0  sec.   Instead,  a  slow  continuous  surface 
reaction  occurs  causing  a  low  regression  rate  during  the 
shut-down  period.   As  a  result,  additional  heat  is  trans- 
ferred to  the  fuel  during  shut-down,  and,  as  compared 
with  the  thermal  profiles  in  Fig.  40a,  profiles  now 
have  less  of  a  slope  and  more  heat  stored  in  the  fuel. 

The  results  of  Case  I  were  considered  an  accurate 
representation  of  the  experimental  firing.   The  major 
purpose  of  the  parametric  study  could  now  be  accomplished: 
to  compare  the  calculated  instantaneous  regression  rate 
behaviors  on  restarts  following  varied  shut-down  durations. 
Three  shut-down  durations  were  considered:   Case  I  with 
a  2.0  sec  shut-down,  Case  II  with  a  0.9  sec  shut-down  and 
Case  III  with  a  4.0  sec  shut-down.   The  instantaneous  re- 
gression rate  behaviors  are  presented  in  two  graphical  forms 

First,  in  Fig.  42,  a  nondimensional  form  similar  to 
that  used  by  Marxman,  Wooldridge  and  Muzzy  in  Ref.  12, 
the  transient  regression  rate  behaviors  for  Cases  I,  II 
and  III  are  compared.   Results  indicate  that  increasing 
shut-down  duration  (up  to  an  undetermined  time)  increases 
the  degree  of  overshoot.   Second,  in  Fig.  43,  if  the 
results  are  presented  in  the  dimensions  of  regression  rate 
and  time,  the  results  are  more  obvious.   The  steady-state 
regression  rates  for  the  radiative  and  nonradiative  systems 
are  shown  in  dashed  and  solid  lines  respectively.   The 
regression  rate  tail-off  associated  with  the  drop  in 
chamber  pressure  should  lead  to  extinction.   However,  slow 
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burning  is  continued  during  the  shut-down.   The  restart 
pressurizations  were  initiated  at  5.9,  7.0  and  9.0  sec 
for  Cases  I,  II,  and  III  respectively.   The  computed  trans- 
ient regression  rates  are  plotted  for  comparisons.   The 
regression  rate  overshoots  are  now  obvious  and  are  clearly 
increasing  with  increased  shut-down  durations.   The  re- 
start of  the  nonradiative  datum  case  is  also  plotted  for 
comparison.   Results  such  as  ignition  delay  and  the  time 
to  reach  90%  of  initial  steady-state  regression  rate  on 
restart  are  tabulated  for  comparison  in  Table  III.   It 
should  be  noted  that  the  time  to  reach  90%  of  the  initial 
regression  rate  on  restart  following  any  shut-down  dura- 
tion is  nearly  one-fifth  of  the  time  to  reach  90%  of  the 
steady-state  regression  rate  on  ignition.   That  is,  steady- 
state  regression  rates  on  reignition  are  achieved  five 
times  as  fast  as  on  initial  ignition. 

VI-D:   Hypothetical  High  Energy  Fuel 

The  physical  properties  of  a  hypothetical  fuel,  whose 
physical  ingredients  were  similar  to  a  fuel  consisting  of 
50%  Boron,  20%  ammonium  perchlorate,  AP ,  and  30%  binder, 
were  substituted  for  the  Plexiglas  properties  for  a  series 
of  calculations  to  determine  the  regression  rate  over- 
shoots that  may  occur  in  a  practical  high  energy  fuel. 
The  fuel  was  assumed  to  be  opaque  to  thermal  radiation. 
Surface  temperature  and  latent  heat  were  varied  for  three 
cases,  and  the  values  for  each  case  is  listed  in  Table  III 
(Case  IV,  V  and  VI) .   The  transient  regression  rates  are 
plotted  in  Fig.  44.   The  regression  rate  overshoots  are 
increased  as  the  latent  heats  are  decreased  and  the  surface 
temperatures  increased. 

An  additional  case,  similar  to  Case  V,  was  considered 
with  a  4.0  sec  shut-down  duration.   As  all  the  other 
initial  parameters  were  equal,  the  regression  rate  transients 
following  ianition  were  similar.   However,  the  regression 
rate  behavior  following  the  2.0  sec  and  4.0  sec  shut-downs 
were  different  and  are  compared  in  Fig.  45.   The  overshoot 
following  the  2.0  sec  shut-down  was  18%  greater  than  steady- 
state,  whereas  the  overshoot  following  a  4.0  sec  shut-down 
was  20%  greater. 
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CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS  WITH  SUGGESTIONS  FOR  FUTURE  WORK 


VII-A:   Summary  and  Conclusions 

From  an  analytical  study  of  the  shut-down  and  restart 
processes,  it  was  concluded  that  if  proper  attention  is 
focused  on  the  inter-relationships  between:   (1)  the  non- 
steady  thermal  state  of  the  fuel  during  shut-down,  (2)  the 
instantaneous  regression  rate  behaviors  following  a  shut- 
down duration  and  (3)  the  pressure  transients  following 
reignition,  the  pressure-time  histories  and  fuel  consump- 
tion following  reignition  may  be  anticipated.   Experimental 
results  have  provided  conclusive  evidence  that  such  re- 
lations between  the  thermal  state  of  the  fuel,  the  instan- 
taneous regression  rate,  and  the  transient  chamber  pressure 
do  exist. 

The  thermocouple  probes,  especially  fabricated  for 
this  study,  provided  accurate  subsurface  temperature  his- 
tories and  consistent  surface  temperature  measurements. 
Thermocouple  studies  showed  that  a  nonsteady,  nonuniform 
thermal  state  exists  in  the  fuel  during  shut-down  and  a 
deeper  thermal  profile  develops.   Chamber  pressures  mea- 
sured during  the  transient  and  chamber  filling  interval 
following  reignition  of  a  preheated  fuel  have  everywhere 
exceeded  their  respective  initial  ignition  transient 
pressures.   From  these  experimental  results,  we  can  infer 
that  the  higher  chamber  pressures  are  a  result  of  the 
higher  rates  of  fuel  mass  addition,  and  therefore  faster 
instantaneous  regression  rates,  caused  by  the  preheating 
of  the  fuel. 

However,  experimental  conditions  have  precluded  any 
measurement  of  the  instantaneous  regression  rates  to 
verify  our  inference.   Therefore,  an  additional  parametric 
study,  using  a  computer  to  simulate  an  experimental  firing 
of  the  02-Plexiglas  hybrid  rocket  motor,  was  performed  to 
determine  the  effects  of  preheating  and  shut-down  on  instan- 
taneous regression  rate  behavior.   The  computer  results 
were  conclusive  and  as  expected;  the  calculated  restart 
regression  rate  transients  had  overshoots  as  large  as  21% 
above  steady-state  values.   The  times  to  reach  90%  of 
equilibrium  chamber  conditions  following  reignition  were 
one-fourth  as  long  as  the  times  on  initial  ignition.   The 
regression  rate  overshoots  increased  as  the  shut-down 
duration  was  increased  for  the  shut-down  intervals  inves- 
tigated (0.9  sec  to  4.0  sec).   The  same  computer  simulation 
was  applied  to  three  hypothetical  high  energy  fuels  to  de- 
termine the  influence  of  the  shut-down  duration  on  reignition 
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regression  rate  transients  in  a  practical  hybrid  fuel. 
The  results  were  similar  to  those  obtained  from  the  PMM 
study. 

The  combined  results  from  the  parametric  computer 
study  and  the  experimental  02~Plexiglas  motor  study  proved 
that  radiant  heating  to  the  subsurface  layers  is  not  negli- 
gible and  cannot  be  ignored  as  previous  researchers  have 
assumed.   Radiant  heating  to  the  subsurface  layers  of  the 
fuel  must  be  considered  in  the  heat  conduction  equation  in 
order  for  experimental  results  to  agree  with  calculated 
predictions.   The  effects  of  radiation  in  the  PMM  system 
increased  the  instantaneous  regression  rates  and  increased 
the  depths  of  the  thermal  profiles  in  the  fuel. 

In  conclusion,  the  effects  of  a  preheated  fuel  (as  a 
consequence  of  a  previous  ignition  and  a  temporary  shut- 
down) will  greatly  influence  the  chamber  filling  interval 
and  the  instantaneous  regression  rates  following  reignition. 
Therefore,  attention  must  be  focused  on  the  overall  thermal 
distribution  during  shut-down  in  order  to  predict  satis- 
factory pressure-time  histories  and  design  the  reliable 
structural  components  capable  of  withstanding  the  excessive 
rates  of  pressurization  following  reignition. 

VII-B:   Suggestions  for  Future  Work 

Future  experimental  work  should  be  directed  at: 

1.  determining  the  upper  limit  of  the  shut-down 
duration   needed  to  eliminate  the  influence  of  the 
preheated  fuel, 

2.  extending  the  computer  simulation  to  predict  the 
time  when  the  preheated  fuel  no  longer  contributes 
to  an  overshoot  on  reignition, 

3.  investigating  the  minimum  time  for  steady  chamber 
conditions  to  be  regained  following  step  changes  in 
the  oxidizer  flow  and  chamber  pressure,  thereby  sim- 
ulating throttling  maneuvers, 

4.  incorporating  the  experimental  results  from  a  simple 
polymeric  fuel  (Plexiglas)  and  the  computer  pre- 
dictions for  a  metal  loaded  fuel  (for  example,  40% 
Al ,  60%  Plexiglas)  to  predict  the  restart  process 
for  a  metal  loaded  hybrid  rocket  motor. 
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Table  II 


Datum  Case  Properties  Used  in  the  Calculations 
for  the  O^-Plexiglas  System 

Thermal  Conductivity,  J$l,    cal/cm  K  sec  0.001213 

Specific  Heat  of  gas,  c  ,  cal/g-c  0.35 

Density,   J>  -,  g/cm  1.683 

Thermal  Diffusivity,  yC     ,  cm  /sec  0.002325 

Ratio  of  Specific  Heats,  "V  ,  1.2 

Average  Molecular  Weight,  combustion  gases  g/g-mole      28.0 

Initial  Temperature,   T-.,   C  2101 

Throat  Diameter  of  Nozzle,   d.,  cm  0.434 

Diameter  of  Motor  Port   d,   cm  2.54 

Maximum  Chamber  Pressure,   P  ,  atmos  24.80 

Surface,  Temperature,   T  ,   C  372.0 

Latent  Heat,   Ah,   cal/g  -351.0 

2 

Igniter  Heat  Flux,  Q,  cal/cm   sec  13065 

Flame  Temperature,   Tf,  C  2732.0 
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Fig.  3   Turbulent  Boundary  Layer  Combustion  Model   (taken 
from  Marxman  and  Gilbert,  Ref.  6) 
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Fig.  5   Steady-State  Thermal  Profile  in  a  Hybrid  Rocket 
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Fig.  6   Transient  Regression  Rate  Behavior  during  the 
Establishment  of  a  Thermal  Profile  in  a  Hybrid 

Fuel   (taken  from  Marxman  et  al.,  Ref.  12) 
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Fig.  22   Log  Tx  vs  Time  for  Run  18  with  Measured  Surface 
Temperatures  Indicated 


Fig.  2  3   Log  T   vs  Time  for  Run  2  2  with  Measured  Surface 
Temperatures  Indicated 
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Fig.  24   Log  T   vs  Time  for  Run  23  with  Measured  Surface 
Temperatures  Indicated 
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Fig.  26   Log  T  vs  Time  for  Run  22  to  Compare  the  Time  to 

Establish  Thermal  Equilibrium,  as  Indicated  by  the 
Straight  Line  Approximation,  with  Time  to  Establish 
90%  of  the  Steady-State  Chamber  Pressure,  t#9Pc 
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FLAME  SPREADING   DURING  HYBRID  IGNITION 
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Fig.    37      Flame   Distribution   During    Steady   Burning   of   Hybrid   Rocket 
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APPENDIX  A 
THERMOCOUPLE  TECHNIQUE 

App.  A-l:   Considerations  of  the  Possible  Errors  in  Thermocouple  Data 

The  1-mil  chromel-alumel  thermocouples  used  in  this 
study  were  purchased  from  Omega  Engineering  Inc. ,  Stamford, 
Connecticut.   These  thermocouples  provided  accurate  re- 
sponses over  the  temperature  range  of  interest  and  small 
enough  wire  diameter  to  enable  temperatures  to  be  meas- 
ured without  greatly  disturbing  the  surrounding  material. 
Yet,  they  were  strong  enough  to  withstand  positioning  on 
Plexiglas  rods  and  embedding  in  the  fuel.   A  Plexiglas 
grain  instrumented  with  four  thermocouples  embedded  in 
the  fuel  is  shown  before  and  after  firing  in  Fig..  46. 

Certain  errors  are  inherent  to  every  temperature 
study.   A  brief  listing  follows  to  describe  some  of  these. 

(1)  The  major  problem  associated  with  any  procedure 
is  the  uncertainty  in  location  of  the  point  at  which  the 
temperature  is  being  measured.   The  uncertainty  arises  from 
two  factors.   The  first  uncertainty  is  locating  the  actual 
position  of  the  thermocouple  at  a  prescribed  depth  from  the 
surface.   Each  successive  measured  dimension  used  to  locate 
the  bead  compounds  the  tolerance  errors  given  on  the  prev- 
ious dimension.   The  second  uncertainty  is  the  error  in  lo- 
cating the  effective  junction  of  the  thermocouple.   The  ef- 
fective junction,  as  defined  by  Brewer,   is  the  point  in 
the  material  whose  temperature  corresponds  to  the  measured 
temperature  output  of  the  thermocouple.   For  example,  the 
recorded  temperature  may  be  assumed  to  be  the  temperature 
on  the  surface  of  the  bead  when  in  fact  it  is  the  tempera- 
tures at  the  center  of  the  bead. 

(2)  Considerable  errors  in  temperature  measurements 
may  result  from  the  disturbance  in  the  material  caused  by 
the  presence  of  the  temperature  probe  itself.   In  particular, 
lead  wires  can  cause  considerable  heat  to  be  conducted  away 
from  the  point  of  measurement.   Studies  of  this  conduction 
error  made  by  Brewer^  show  that  the  errors  are  significant 
if  the  thermocouple  lead  wires  near  the  bead  are  parallel  to 
the  direction  of  heat  flow.   Strittmater,  et  al.°  agree 
that  it  is  attractive  to  consider  directing  the  leads  along 
an  isotherm  to  prevent  heat  losses  through  them.   However, 

as  the  authors^  indicate,  an  isotherm  could  be  difficult 
to  locate  due  to  nonuniform  burning  surfaces  resulting  from 
propellant  roughness.  Also,  as  the  entire  surface  is  seldom 
ignited  all  at  once,  placement  of  the  lead  wires  parallel  to 
the  initial  grain  surface  involves  the  risk  that  the  burning 
surface  will  reach  the  lead  wire  before  the  bead  and  preheat 
it  by  conduction. 

A  compromise  between  the  two  considerations  mentioned 
above  was  made  for  this  study.   Lead  wires  in  the  vicinity 
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of  the  bead  were  oriented  perpendicular  to  the  heat  flow  to 
minimize  the  conduction  error.   At  a  sufficient  distance 
away  from  the  bead,  the  lead  wires  were  directed  away  from 
the  surface  plane  to  avoid  their  heating  before  the  bead. 

(3)   The  temperature  recorded  by  a  sensing  probe  is 
dependent  on  its  surrounding  environment.   The  influence 
of  radiative  heating  can  produce  large  increases  in  temp- 
erature responses  because  the  radiation  effect  is  propor- 
tional to  the  temperature  of  the  source  raised  to  the  fourth 
power  minus  the  temperature  of  the  sink  raised  to  the  fourth 
power.   However,  it  was  initially  assumed  for  these  experi- 
ments that  the  radiation  effect  on  the  thermocouple  outputs 
would  be  insignificant  and  could  be  ignored.   This  assump- 
tion was  based  on  two  factors:   1)  the  emissivities  of  PM 
in  the  temperature  range  of  interest  is  close  to  unity1^ 
and  2)  the  temperature  sensing  probes  were  made  from  mater- 
ials that  would  reflect  most  of  the  radiant  energy  reaching 
the  bead.   Therefore,  as  previous  workers  have  assumed^' 15, 18^ 
the  radiation  effects  would  be  negligible.   However,  as  dis- 
covered from  the  results  of  the  thermocouple  data,  radiation 
to  the  subsurface  regions  of  the  fuel  is  significant.   Thus, 
a  solution  to  the  heat  conduction  equation  considering  only 
convective  heating  would  yield  poor  results  when  compared  to 
experimental  data,  but  when  two  radiation  parameters  are  con- 
sidered, the  solution  produced  favorable  results.   Therefore, 
it  has  be£n  an  incorrect  assumption  by  previous  researchers 
to  disregard  the  radiative  influence  on  temperature  measure- 
ments in  Plexiglas  systems. 
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App.  A-2:   Thermocouple  Mounting  Procedure 

A  mounting  procedure  was  devised  to  minimize  all  the 
errors.   The  thermocouples  were  positioned  and  maneuvered 
with  the  aid  of  a  micromanipulator  and  a  Stereo-zoom  micro- 
scope.  The  procedure  used  for  mounting  and  installing  are 
discussed  in  the  following. 

(1)  The  thermocouple  was  placed  on  a  1/4  inch  diam- 
eter rod  in  a  machined  groove  0.005  inches  in  depth  across 
the  diameter  of  the  face.   Plexiglas  rods  1/4  x  2  inches 
were  used  so  as  not  to  effect  the  properties  of  the  fuel. 
The  thermocouple  was  positioned  in  the  groove  with  the  bead 
as  close  to  center  as  possible.   A  drop  of  PS-18,  a  plexi- 
glas cement,  was  used  to  fasten  the  thermocouple  bead  to 
the  rod.   This  cement  was  used  since  it  hardens  to  the  same 
material  as  the  surrounding  fuel.   Therefore,  it  should  not 
effect  ablation  in  the  areas  of  the  bead  or  the  fuel  sur- 
rounding the  bead. 

(2)  The  lead  wires  were  also  cemented  across  the  face 
of  the  1/4  inch  rod.   When  embedded  in  the  fuel,  the  lead 
wires  will  be  perpendicular  to  the  heat  flow  in  the  vicinity 
of  the  bead  and  thereby  minimize  the  conduction  error  de- 
scribed earlier.   The  remaining  lengths  of  lead  wire  were 
drawn  along  the  sides  of  the  rod  and  cemented  into  0.132 
inch  groove.   An  Omega,  1/4  inch  Chromel-Alumel,  adapter 
was  fitted  to  the  end  of  the  rod  and  the  lead  wires  were 
connected  to  their  respective  terminals.   The  rod,  complete 
with  thermocouple,  was  ready  for  installation  in  the  fuel. 

(3)  Radial  holes  with  flat  bottoms  to  assure  positive 
depth  locations  were  drilled  into  the  Plexiglas  fuel  grains. 
Normally  each  hole  was  measured  to  locate  each  thermocouple 
at  a  different  depth.   The  holes  were  half  filled  with  the 
same  Plexiglas  cement,  PS-18,  and  the  rods  were  carefully 
inserted,  forcing  some  of  the  cement  and  most  of  the  air 
bubbles  to  the  surface.   This  procedure  assured  a  gas-tight 
seal  near  the  bead.   The  rods  were  secured  with  tape  and  at 
least  24  hours  was  allowed  for  drying  before  firing  the  in- 
strumented grain. 

The  mounting  procedure,  although  an  arduous  task  re- 
quiring enduring  patience,  produced  a  reliable  thermocouple 
probe.   One  hundred  percent  of  the  thermocouples  installed 
in  the  fuel  functioned  property. 
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APPENDIX  C 
OPERATING  PROCEDURE 
FOR 
HYBRID  MOTOR 
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A.  Pre  Firing  Procedure 

1.  Turn  on  all  electrical  equipment  to  allow  for  30  min.  warm  up. 

2.  Select  and  install  0_  orifice. 

3.  Select  and  install  CH.  orifice. 

4.  Weigh  and  install  graphite  insert. 

5.  Install  pressure  transducers:   include  RTV  \   disk  where  needed, 

6.  Clean,  weigh,  and  install  grain. 

7.  Install  nozzle  plate. 

8.  Connect  all  recording  devices. 

9.  For  Electrical  hookup  procedure  go  to  Section  B. 

10.  Turn  on  main  N«  pressure  (Mechanical  Room) • 

11.  Open  N_  control  pressure  valve  (Test  Cell)  . 

12.  For  calibration  procedure  go  to  Section  C. 

13.  Open  CH.  bottle  valve  and  set  CH.  regulator. 

14.  Open  0«  bottle  valve  and  record  pressure. 

15.  Plug  in  ignitor  control  box  at  wall  outlet. 

16.  Plug  in  spacing  grid  hamp. 

17.  Ignitor  control  safe  switch  to  on. 

18.  Install  exhaust  duct. 

19.  Remove  keys  and  clear  area  #Z. 

20.  Open  blower  vent  and  vacate  test  cell. 

21.  Turn  on  blower. 
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B.  Electrical  Procedure 

1.  Plugs  from  test  area  #Z  (4  pins)  to  transducer  input  #1. 

2.  Calibrate  voltmeter  @  9000. 

3.  Zero  voltmeter  @  0000. 

4.  Zero  amplifiers  @  10  volts  with  shunts. 

5.  Calibrate  transducers  with  10  volts  thru  input  conductioners 
using  volt  button. 

6.  Balance  transducer  thru  amplifier  output  with  balance  buttons 
to  zero. 

7.  Connect  marker  blips  using  reversed  (bronze)  cords. 

8.  Fl  to  CH2. 

9.  01  to  CH4. 

C.  N2  Calibration 

1.  From  A  #11. 

2.  Check  N2  line  supply  from  area  4. 

3.  Turn  on  24  volt  power. 

4.  N2  supply  switch  on  panel  to  on. 

5.  Calibrate  @  pressures. 

6.  Dump  N2  pressure  from  panel  switch. 
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D.  Firing  Procedure 

1.  Turn  on  circuit  breaker  14,  17,  12,  and  DC  supply. 

2.  Turn  on  N_  supply  switch  on  panel. 

3.  Set  Op  reservoir  pressures. 

4„  Check  ignitor  mode  switch,  ignitor  timer  and  run  duration  timer, 

5.  Reset  automatic  0?  switch  to  down  position. 

6.  Power  to  area. 

7.  Warning  lights  on. 

8.  Siren  on. 

9.  Firing  mode  switch  to  READY. 

10.  Firing  sequence. 

11.  Terminate  sequence  with  reset  0«  switch  to  up  position. 

E.  Shut  Down  Procedure 

1.  Np  purge  for  2  sec. 

2.  Siren  off. 

3.  Firing  mode  switch  to  off. 

4.  Remove  keys  from  console. 

5.  Turn  off  circuit  breaker  #17. 

6.  Blower  on  for  15  sec. 

7.  Enter  test  cell  and  turn  off  ignitor  safety  switch. 

8.  Remove  ignitor  control  box  plug  from  wall  outlet. 

9.  Close  CH.  and  0«  bottle  valves. 

10.  Vacate  test  cell  and  turn  blower  on. 

11.  Insert  keys. 

12.  Firing  mode  switch  to  SAFE. 
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13.  Bleed  02  line  and  reservoir  pressure. 

14.  Bleed  CH.    line  intermittent  N2  purge. 

15.  Firing  mode  switch  to  off. 

16.  Remove  keys  from  console. 

17.  Turn  off  N2  control  pressure  valve. 

18.  Insert  keys. 

19.  Firing  mode  safety  switch  to  safe. 

20.  Bleed  purge  pressure. 

21.  Console  N2  supply  switch  to  off. 

22.  Warning  lights  off. 

23.  24  volt  DC  power  off. 

24.  Circuit  BK  #14  off. 

25.  Break  down  and  clean  motor. 

26.  Clean  up  t&st  cell. 
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